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Within the non-perturbative 1/N expansion, we discuss numerical methods for calculating multi-loop Feynman 
graph needed to derive physical scattering amplitudes. We apply higher order 1/N methods to the scalar sector 
of the standard model, and show the existence of a mass saturation effect. The mass saturation has direct 
implications for future searches at the LHC and at possible muon colliders. 



1. INTRODUCTION 

There has been a lot of work lately in extending 
the calculation of radiative corrections to higher 
loop-orders and to multi-leg processes. The pro- 
ceedings of this Workshop are suggestive for the 
complexity this field has reached since the pio- 
neering work of Martinus Veltman and Gerardus 
't Hooft. Higher order QCD radiative corrections 
are necessary for interpreting the experimental 
data which the LHC experiments will provide. 
Higher order electroweak corrections are needed 
for interpreting precision electroweak measure- 
ments. Due to the time frame for the construc- 
tion of future colliders, precision measurements 
will continue to be a main tool for new physics 
searches. 

In this contribution we would like to discuss the 
use of multi-loop Feynman graphs for recovering 
results which are normally not within the scope of 
perturbation theory, namely deriving scattering 
amplitudes in the case of strongly coupled field 
theories. 

A central question in particle physics is how 
is the electroweak symmetry broken. The ex- 
perimental data currently available is compatible 
with a standard model scalar sector, and favors a 
light Higgs boson. However, additional degrees of 
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freedom beyond the standard model have the po- 
tential to shift these fits. It is not at all excluded 
that the LHC experiments will see a scalar reso- 
nance at a considerably higher energy. 

At strong coupling, the scalar sector was in- 
sufficiently explored. Qualitatively, one may ex- 
pect new phenomena to appear: a Higgs parti- 
cle strongly coupled to the vector bosons and to 
itself, anomalous vector boson self-interactions, 
and possibly the appearance of a spectrum of ad- 
ditional resonances in the scalar sector. Pertur- 
bative calculations loose their predictive power 
when radiative corrections blow up in higher loop 
orders, and the dependency on the renormaliza- 
tion scheme becomes substantial [ . 

Ideally, one needs a solution for the scalar sec- 
tor which is valid at strong coupling as well as at 
weak, and which is free of renormalization scheme 
uncertainty. This can be accomplished by explic- 
itly summing up all loop orders within a non- 
perturbative 1/N expansion. If a solution of suf- 
ficient accuracy is desired, such that it can be 
used in phenomenological calculations to compete 
with ordinary perturbation theory, then the 1/N 
expansion must be carried out at higher order. 

2. THE 1/N SOLUTION 

We are interested in a solution of the scalar sec- 
tor of the standard model when the self-coupling 
of the Higgs field becomes strong. The non- 
perturbative effects are thus given by the scalar 
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field self-interaction, while the gauge coupling re- 
mains perturbative. It is then natural to re- 
sort to the equivalence theorem to relate am- 
plitudes involving longitudinal electroweak gauge 
bosons to the corresponding amplitudes involving 
would-be Goldstone bosons. The problem is be- 
ing reduced to solving a linear sigma model non- 
perturbatively. 

2.1. The auxiliary field formalism 

It is sufficient to consider the scalar sector 
alone, which is an SU (2)-symmetric linear sigma 
model. In the following, the scalar sector is ex- 
tended to an 0(iV)-symmetric sigma model - 
with the standard model case recovered for N = 4 
— such that an expansion in powers of 1/N can 
be performed: 
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From the Lagrangian above, it is easy to derive 
scattering amplitudes at leading order in 1/N. 
Typically, this involves summing up a geomet- 
ric series of one-loop bubble diagrams which have 
Goldstone bosons in the loop. 

However, the auxiliary field formalism which 
was proposed in ref. [ ^| proves to be most useful 
in organizing sub- leading orders in the 1/N ex- 
pansion in a diagrammatically manageable way. 
It consists in introducing an additional unphysi- 
cal field x m the Lagrangian: 



Co = 



3N \q 2 22 
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Xo + const. 
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Because the equation of motion for the auxil- 
iary field x is a constant, this addition does not 
change the dynamics. Green's functions having 
Higgs and Goldstone bosons on the legs are the 
same, whether calculated with the Feynman rules 
given by L\ or by £2- 



The Feynman rules, on the other hand, are 
changed. C\ contains the trilinear and quartic 
vertices <77T7r, aaa, aaaa, acririr, and tttttttt. Here 
a and ir are the massive and massless modes stem- 
ming from Lagrangian C\ after spontaneous sym- 
metry breaking, respectively. £2 contains only 
the trilinear couplings x a<J an d X 7r7r - 

This simplifies enormously the topological clas- 
sification of Feynman diagrams according to their 
power of 1/N beyond leading order. The reader 
can easily convince himself of the utility of the 
auxiliary field formalism by writing down the di- 
agrams which contribute to Goldstone-Goldstone 
scattering at NLO in 1/N in both formalisms. 

2.2. Tachyonic regularization 

By summing up the chains of one-loop bubble 
self-energy insertions, one encounters an ultravi- 
olet renormalon. This gives rise to an additional, 
tachyonic pole in the propagators, apart from the 
expected physical spectrum containing one Higgs 
boson and N — 1 Goldstone modes. 

By direct evaluation of the leading order con- 
tribution in 1/N to the two-point functions, one 
obtains the following propagators [ H, E| : 
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D xx {s) 
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where 
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Here a^(s) is the ultraviolet finite part of 
the one-loop self-energy bubble diagram, with a 
Goldstone boson in the loop, fj, is the ultraviolet 
subtraction scale. 
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The propagators contain an Euclidian pole in 
the ultraviolet region at an energy s 
by the following transcendental equation 



given 



1 



A? 



32tt : 



log 
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The tachyon scale is in the ultraviolet region for 
low values of the coupling, and tends to move to- 
wards low energy when the coupling is increased. 

In higher order 1/N corrections, the Euclid- 
ian pole appears in the loop momentum integra- 
tion. This induces causality violating contribu- 
tions, even though these effects are numerically 
small as long as the tachyon scale is high enough. 
For this reason, it is necessary to investigate the 
origin of this singularity and find a way for deal- 
ing with it. 

At any finite order of perturbation theory, the 
tachyon pole does not exist. It is an artifact 
of the bubble diagram summation. In the pro- 
cess of summing up all loop orders, an ambigu- 
ity is present, however. The summation of the 
perturbative series determines the result only up 
to functions which vanish in perturbation the- 
ory, such as e -1 / A . Because the residuum of the 
tachyonic pole is precisely such a function which 
vanishes in perturbation theory, its presence or 
absence is completely arbitrary and cannot be 
determined within perturbation theory. For this 
reason, it is justified to restore causality by min- 
imally subtracting the tachyon at its pole since 
the original information stemming from Feynman 
diagrams remains unchanged [ Ef : 



D xx (s) 
D xa {s) 



1 



s - 

IS 

Nv 2 
i 



Nv 



m 2 (s) s- 

m 2 (s) 
s — m 2 (s) 
m 2 (s) 
— m 2 (s) 



•A? 



+ A 2 
K A 2 



A? 



(0) 



were k = [1 + A 2 / (32n 2 v 2 )]^ 1 is the residuum of 
the tachyonic pole. 

The tachyonic regularization can be seen as a 
prescription for summing up the perturbative se- 
ries in a way which preserves causality. 



2.3. 1/N renormalization 

Performing renormalization at NLO in the 1/N 
expansion involves the treatment of ultraviolet di- 
vergences of diagrams with various numbers of 
loops. One way of keeping track of various loop 
order counterterms is to group them into 1/N 
counterterms. Each of the 1/N order countert- 
erms SX, St, 6t x , 5Zn,a, x is a power series in the 
coupling constant A [ Ej : 
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In principle, these counterterms are sufficient 
to absorb the UV divergences of all diagrams of 
NLO in 1/N. 

When calculating 1/N diagrams at NLO, we 
resort to numerical integration. The ultraviolet 
divergence of the diagram must be removed be- 
fore the numerical integration can be carried out, 
but the explicit e expansion and isolation of UV 
poles is cumbersome because of the complexity 
of the diagrams. For this reason we perform a 
BPHZ-type renormalization [ ||] . We subtract the 
divergences of the diagrams according to the for- 
est formula, as shown in figure 1 . The subtracted 
expressions, being finite in the ultraviolet, can be 
integrated numerically. For a given physical pro- 
cess, the renormalization program then means to 
combine all UV subtraction terms from various 
diagrams with each other. Most subtractions can- 
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Figure 1. Ultraviolet subtractions of all NLO 
1/N two- and three-point diagrams of the sigma 
model. The blob on the propagators denotes the 
resummed chain of one-loop Goldstone bubble di- 
agrams. The solid line denotes Goldstone fields, 
the wavy line is the Higgs field, and the dashed 
line is the auxiliary field. 



eel with each other, and the small remaining set 
of subtractions is trivially absorbed into the 1/N 
local counterterms above. We have checked that 
all UV subtractions are polynomial, as they ought 
to be in order to be absorbed into local countert- 



2.4. Numerical solution 

The finite expressions depicted in figure 1 arc 
calculated by numerical integration because so far 
no analytical approach exists for dealing with this 
type of diagrams. 

For a numerical solution, it is advantageous to 
identify first all loop integrations which can be 
performed analytically. These are associated with 
closed Goldstone loops. For the diagrams shown 
in figure 1, it is possible to perform analytically all 
integrations except for one final loop integration. 
The final integration involves the resummed and 
tachyonically subtracted propagators of eqs. 6 
and form factors from one-loop triangle or box 
diagrams involving massless Goldstone bosons in 
the loop. The final loop integration needs to be 
performed numerically. It can be reduced to a 
two-fold integral with the methods of refs. [ |[ H] , 
which do not resort to Feynman parameters. 

2.5. Scheme independence of physical pre- 
dictions 

For a given order in the expansion parameter 
1/N, the Feynman diagrams of all loop orders 
are explicitly summed up. For this reason, the 
final result for a physical scattering amplitude is 
free of renormalization scheme dependence. In 
usual perturbative calculations, a renormaliza- 
tion scheme ambiguity is present because of the 
truncation of the perturbative expansion; this 
ambiguity is of higher order in the coupling con- 
stant. 

We note that subtracted diagrams of the type 
shown in figure 1 are not individually free of 
renormalization scheme dependency. The inter- 
nal subtractions which are performed for making 
the integrals finite in the ultraviolet, are calcu- 
lated at a given subtraction point. This defines 
an intermediary renormalization scheme. 

This renormalization scheme dependence can- 
cels out only in physical results such as scatter- 
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ing amplitudes. In the following section we give 
the results for two scattering processes. We have 
checked explicitly that the two functions involved, 
fx and ji of eqns. 9, are independent of the sub- 
traction point. 

3. COLLIDERS AND THE HIGGS 

MASS SATURATION EFFECT 

Recently, feasibility studies for muon collid- 
ers attracted quite some attention. A muon col- 
lider would be an ideal s-channel Higgs factory. 
For a heavy Higgs boson, two processes domi- 
nate: u + u~ —tH^tt and u + u~ —tH^ 



W£W£,Z L Z L . 

Within the 1/N expansion, the amplitudes for 
these processes at NLO are given by the following 
expressions [ [!): 

1 



M ff 



M 



zz 



s — m 2 (s) [l — 
m 2 (s) 1 



^A(*)] 



(8) 



Here, the correction functions f\ and fi are given 
by a combination of the two- and three-point 
functions defined in figure 1 (see [ ) : 
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In figure 2 we give numerical results for these 
expressions. We plot the shape of the Higgs reso- 
nance for various strengths of the coupling. From 
the positions of the resonance maxima, it can be 
seen clearly that when the coupling is increased, 
a mass saturation effect appears. The position 
of the resonance's peak does not increase above 
a saturation value just under 1 TeV. The pre- 
cise resonance shape and implicitly the saturation 
value are process dependent. 

For comparison, we indicate in figure 2 the 
corresponding peak maxima which would be ob- 
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Figure 2. The Higgs line shape at a Higgs 
factory. We marked the position of the maxima 
of the resonances (solid line for the 1/N result 
and dashed line for the perturbative result at two- 
loop). 



tained by using usual perturbation theory at 
NNLO. A saturation is present qualitatively in 
this curve, too. However, it should be noticed 
that the perturbative curve is affected by large 
radiative corrections and large scheme uncertain- 
ties at such large coupling, and therefore is not 
reliable in the saturation region. 

4. SATURATION EFFECT AT THE LHC 

The main Higgs production mechanism at the 
LHC is gluon fusion which proceeds through a 
top loop [ 0| . The existing perturbative analy- 
ses indicate that the vector boson fusion becomes 
competitive at an energy of the order of 1 TeV. 

In the presence of the saturation effect, the ob- 
servation of a heavy Higgs resonance at the LHC 
will be different. In the saturation zone, the mass 
of the resonance remains more or less the same 
while the width increases with the coupling. The 
resonance becomes flatter and more difficult to 
detect. 

In figure 3 we show a plot of the Higgs width 
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Figure 3. The current knowledge of the Higgs 
width at strong coupling from perturbation the- 
ory an from the 1/N expansion. The parame- 
ters Mpeak and Tpeak are extracted from the 
position and the height of the Higgs resonance 
in fermion scattering as if the resonance was of 
Breit- Wigner type. For the perturbation theory 
curves we give the corresponding values of the on- 
shell mass parameter to#. 



as a function of the Higgs mass. The width and 
mass used in these plots, Mpeak and Tpeak, 
are defined from the line shape of the resonance 
as seen in fermion-fermion scattering. They are 
derived from the position and height of the reso- 
nance as if it were of Breit- Wigner type. In this 
picture we show the calculations available so far 
in usual perturbation theory (LO, NLO, NNLO) [ 
§, and in the 1/N expansion (LO and NLO) [[|. 
It can be seen that the two expansions converge 
nicely towards each other and display a mass sat- 
uration. 

We studied the discovery potential of the LHC 
in the presence of the saturation effect. In ref. 
[ |^] we performed a Monte-Carlo simulation for 
the LHC. Only leptonic channels were considered, 
namely the "golden plated" channel (l + l~)(l + l~), 
and the l + l~vv channel. When neutrinos in the 
final state are involved, the Higgs resonance ap- 
pears as a Jacobian peak in the distribution of 



missing p?. We included the gluon fusion pro- 
cess together with the relevant background. The 
strong interacting Higgs correction was included 
by using the NNLO perturbative calculation - it 
results into a faster Monte-Carlo while being a 
fair approximation of the non-perturbative 1/N 
result. 

We used usual assumptions about the LHC en- 
ergy and luminosity, and asked for a 5a effect 
with respect to the background. The discovery 
potential estimated for the "golden plated" chan- 
nel corresponds then to an on-shell Higgs mass of 
830 GeV. The missing pt channel reaches up to 
an on-shell Higgs mass of 1030 GeV. Note that 
these values for the on-shell mass deviate consid- 
erably from the actual position of the resonance. 
The actual peak can be read from figure 3, where 
the on-shell mass is mapped onto the saturation 
curve. 

5. CONCLUSIONS 

The combinatorial structure of the sigma model 
makes possible an explicit calculation of the non- 
perturbative 1/N expansion at NLO. 

The 1/N solution is valid at strong coupling 
as well as at weak coupling. It is also indepen- 
dent of the intermediate renormalization scheme 
which is being used. At NLO it also provides 
for ultraviolet finite wave function renormaliza- 
tion constants. 

The non-perturbative solution of the sigma 
model is relevant for the standard model via the 
equivalence theorem. It implies that a mass satu- 
ration effect is present in the scalar sector. When 
the Higgs self-interaction is increased, the mass 
of the resonance increases only up to a maximum 
value under 1 TeV, while the width increases con- 
tinuously. 

An interesting question is how cutoff effects 
stemming from new physics at higher energies can 
modify the 1/N solution. This clearly deserves 
further investigation. 
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